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Abstract 

The relationship between MS and pole masses of the vector bosons Z and W is 
calculated at the two-loop level in the Standard Model. We only consider the purely 
bosonic contributions which represent a gauge invariant subclass of diagrams. All cal- 
culations were performed in the linear gauge with three arbitrary gauge parameters 
utilizing the method of asymptotic expansions. The results are presented in analytic 
form as series in the small parameters sin 2 0w arid the mass ratio m 2 z /m 2 H . We also 
present the corresponding on-shell mass counter-terms for the massive gauge bosons, 
which will be needed for the calculation of observables at two-loops in the on-shell 
renormalization scheme. 



1 Introduction 



Precision Physics of the electroweak gauge bosons Z and W started about 12 years ago at 
the LEP storage ring with the ALEPH, DELPHI, L3 and OPAL experiments and ended 
just recently with the dismantling of the LEP installation. In particular the very accurate 
determination of the masses and the couplings to the fermions revealed unexpectedly rich 
information about the quantum correction of the Standard Model (SM) [0. Calculations 
of higher order corrections thus gained increasing importance. At the one-loop level these 
calculations for the relevant "2 fermions into 2 fermions" processes were completed before 
LEP started operating in 1989 B. These SM predictions enabled an indirect determination 
of the top mass which culminated in the top discovery at the Tevatron. Now after the top 
mass has been fixed with rather good accuracy, the indirect bound to the Higgs mass, the 
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last missing SM parameter, is the main goal. The knowledge of the actual value of the Higgs 
mass is extremely important because it determines how Higgs physics will look like at future 
colliders like the LHC or TESLA [[J. Since the sensitivity of SM predictions on the Higgs 
mass is weak the precise meaning of the indirect Higgs mass bounds depend crucially on the 
accuracy of the theoretical predictions. Fortunately a lot of important theoretical progress 
has been made in the last decade with the calculation of leading and some sub-leading two- 
loop effects [||-@. However, no complete two-loop calculation could be achieved so far, 
because such calculations are hampered by the dramatic increase in complexity encountered 
in such calculations. How important the precise evaluation of radiative corrections is may 
be illustrated by the following fact: taking only the leading corrections, the shift Aa em in 
the fine structure constant and the quadratic top mass correction Ap top in the relationship 
between neutral and charged current effective couplings, predictions are about 10a off from 
the data for most of the precisely known observables like sin 2 #f ff or Mw [10 , |TT| • Thus 
the sub-leading effects are huge in relation to current experimental precision. Therefore the 
issue of sub-leading two-loop corrections has to be taken very seriously. They easily may 
obscure the interpretation of the indirect Higgs mass bound obtained from LEP experiments 
by using SM predictions which are incomplete at the two-loop level. 

Although we are a long way from filling the g first step we calculate the full 

bosonic two-loop electroweak corrections to the on-shell masses of gauge bosons. In other 
words, we evaluate the two-loop renormalization constants, the relations between bare, MS 
and on-shell masses, which are required as an essential part for the two-loop renormalization 
program in the on-shell scheme, with the fine structure constant a and the gauge boson 
masses M w and M z as independent input parameters. 

The programs developed for performing the present calculations can be applied without 
modifications for the calculation of the on-shell wave-function renormalization constants of 
the gauge bosons. Since these quantities are gauge dependent we don't reproduce corre- 
sponding results here. As yet we have not calculated any observable quantity. However, the 
new set of corrections will be important for future, complete electroweak two-loop calcula- 
tions of physical observables. Also after the shutdown of LEP it is important to continue such 
calculations because the question how additional corrections affect the Higgs mass bound 
can be answered retrospectively, once given the precise LEP/SLC results. Full two-loop 
calculations would be indispensable in future in any case if a project like TESLA with the 
GigaZ option would be realized. 

The most recent essential progress here was achieved in the calculation of the top-quark 
contributions to the two-loop electroweak corrections. The corresponding contributions to 
the p-parameter were considered in ||, the one's to Ar, which determines the Mw — Mz 
relationship, in |7j and to the Z boson partial widths in Two different approaches — 
the asymptotic expansion method [12|| and numerical integration) 13 1 — have been used to 



perform these calculations. One of the important steps when performing these calculations 
is the two-loop renormalization of the gauge boson masses, which also is contributing to the 
sin 2 9w renormalization [H-[|TH]. In the SM so far no complete analytical calculation of the 



5 Another important step forward was the calculation of the two-loop QED correction to the muon decay 
width [61. 
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two- loop renormalized propagator has been carried out [fT7|] . The first available results were 
given for zero external momentum ||, when the original diagrams may be reduced to a set of 



bubble- type integrals with different masses [18]. In |19[ the two-loop unrenormalized fermion 



corrections to the gauge boson propagator have been presented for off-shell momentum in 
the general linear gauge. The results are presented in terms of scalar master integrals 
with several different mass scales, the masses of fermions and bosons. For the evaluation 
of these master diagrams analytical results |2(| ^] or one-fold integral representations are 



available 22 



The main aim of the present paper is to present a calculation within the Standard Model 
of the two-loop bosonic contributions to the relationship between MS and on-shell masses 
of the gauge bosons (W, Z). This relationship, alternatively, will be represented in terms of 
on-shell renormalization constants. We shall discuss in some detail the algorithm used for 
calculating two-loop electroweak corrections for on-shell quantities in an arbitrary gauge. 

The paper is organized as follows. In Section 2 we briefly reconsider the definition of 
the pole mass of the massive gauge bosons within the Standard Model. The calculations 
have been performed with the help of computer programs which will be described in some 
detail in Section 3. In Section 4 we discuss the UV renormalization of the pole mass and the 
interrelation of our results with the standard renormalization group approach. In particular 
we make several cross-checks of the singular 1/e 2 - and 1/e-terms. The numerical results 
for the finite parts are presented and discussed in Section 5. For further technical details 
and some useful formulae we refer to four appendices. In Appendix A we collect results 
for the one-loop propagator type diagrams. Special attention is given here to the e-parts 
of the corresponding integrals which are needed for the two-loop calculation. Appendix B 
and C collect one-loop results in d = 4 and d ^ 4, respectively. They are included for 
completeness. In Appendix D we present the analytical coefficients which are the main 
results of our investigation. 



2 Pole mass 

The position of the pole sp of the propagator of a massive gauge boson in a quantum field 
theory is a solution for p 2 at which the inverse of the connected full propagator equals zero, 
i.e., 

sp-m 2 -U(sp,m 2 ,---) =0, (2.1) 

where H(p 2 , • ■ ■) is the transversal part of the one-particle irreducible self-energy. The latter 
depends on all SM parameters but, in order to the keep notation simple, we have indicated 
explicitly only the dependence on the external momentum p and in some cases also m, where 
m is the mass of the particle under consideration. This can be either the bare mass m or 
the renormalized mass defined in some particular renormalization scheme. 

Generally, the pole sp is located in the complex plane of p 2 and has a real and an 
imaginary part. We write 

s P = M 2 - iMV. (2.2) 
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The real part of ( |2.2| ) defines M which we call the pole mass^, while the imaginary part is 
related to the width T of the particle. This is the natural generalization of the physical mass 
of a stable particle, which is defined by the mass of its asymptotic scattering state. 

For the remainder of the paper we will adopt the following notation: capital M always 
denotes the pole mass; lower case m stands for the renormalized mass in the MS scheme, 
while m denotes the bare mass. In addition we use e and g to denote the U(l) em and 
SU(2) L couplings of the SM in the MS scheme. 

In perturbation theory (|2.1|) is to be solved order by order. To two loops we have the 



following solution of (|2.1|) 

Sp = m 2 + n (1) (m^m 2 ,---)+n (2) (m 2 ,m 2 ,---) + n (1) (m 2 ,m 2 ,---)n (1), (m 2 ,m 2 ,---), (2.3) 

which yields the pole mass M 2 and the width T at this order, LT^ is the bare (m = mo) or 
MS -renormalized [m the MS -mass) L-loop contribution to II, and the prime denotes the 
derivative with respect to p 2 . In this way we need to evaluate propagator type diagrams and 
their derivatives at p 2 = m 2 . Diagrammatically the self-energy contributions are shown in 
Fig.0. 

The (on-shell) mass counter terms SM 2 for the on-shell renormalization scheme are ob- 
tained by considering the real part of ( |2.3| ) upon identifying the r.h.s with the bare expression, 
setting m 2 = ml = M 2 + SM 2 and solving for SM 2 = (SM 2 )^ + (5M 2 )^ to two loops. 

In this paper we show by explicit calculation at the two-loop level that the bosonic 
correction to the pole sp (and hence to the on-shell mass counter-term) is a gauge invariant 
and infrared stable quantityQ The free propagator of a massive vector boson in the linear 
i?g gauge reads 



i 



D M = -ri—i | -9r + (1 - OzjT^L* ) > ( 2 - 4 ) 



p2 _ m 2 I pi _ ^ m 2 



where £ is a gauge parameter. We decompose the vector boson self-energy U lJiU (p 2 ,m 2 , 
into a transverse H{p 2 ) and a longitudinal L(p 2 ) part 

U, u (p 2 ,m 2 , ■■■) = (g» v - ^) n(p 2 ) + P -^L{p 2 ). (2.5) 

In ( |2.3|) only the transverse part contributes and the dressed propagator reads 

^w- p2 _ m 2_n( P 2 ) p 2 ) + p 2 "' { > 

The simple relation between the full propagator and the irreducible self-energy only holds if 
there is no mixing, like for the iy-boson. In the neutral sector, because of 7 — Z mixing, we 



6 Throughout this paper we identified the terms pole mass and on-shell mass. 

7 From general considerations we know (see [|l6], for example), that after summing all diagrams, including 
the tadpoles, at a given order of perturbation theory the location of the pole of the propagator is a gauge 
invariant quantity. The gauge invariance of the 2-loop massless fermion correction to the pole mass of the 
gauge bosons was verified in |19| . 
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cannot consider the Z and 7 propagators separately. They form a 2 x 2 matrix propagator, 
so that 02.11) is modified into (see details in Hl6|, 1231) 



s P - ml - U zz (sp) ~ \ {SP s = 0. (2.7) 

s P - ll 77 (s P ) 

We note that the U 2 Z mixing term starts to contribute at the two-loop level. Obviously, 
we do not need to compute II 77 here since it starts to play a role only beyond the two-loop 
approximation. In the sequel we will denote the self-energies by HV (V = W, Z) with 



and 



Rw(p 2 , ■■■) = Rww(p 2 



p 2 - n 77 (p 2 , • ■ ■) ' 

Thus, formally, the form (|2.1|) applies for both the W and the Z. 

The non-zero imaginary part (width) ( |2.2p of the on-shell gauge boson self-energy appears 
as soon as the fermions are included. For the bosonic contributions alone the imaginary part 
of n(p 2 ) on the mass-shell is zero at the two-loop level (see below). 



3 Program part 

In order to find the relations between the pole masses Mf , and the MS masses m 2 z , m^y 
we have to compute the one- and two-loop self-energies for Z- and P^-bosons at p 2 = m 2 z 
and p 2 = miy, respectively. The complete set of topologies that occurs in this calculation is 
shown in Fig. [I]. In order to be able to work with manifestly gauge parameter independent 
renormalization constants we have to include the Higgs tadpole diagrams. 

While at one-loop order we have about 50 diagrams, in the two-loop approximation the 
number of diagrams is about 1000, which requires an automatized generation and evaluation 



of diagrams. We use QGRAF |23] to generate the diagrams and then the C-program 
DIANA [BIJ to produce for each diagram an input suitable for our FORM |26| packages^. 
For two-loop propagator type diagrams with several masses a complete set of recurrence 



relations is given in p7| . It allows us to reduce all tensor integrals to a small set of so-called 



master-integrals. However, the master-integrals which show up in the SM are not expressible 
in terms of known functions but may be written e.g. as one- fold integrals ||22| . Instead of 
using these explicit formulae we resort to some approximations here, namely, we perform an 
appropriate series expansion in mass ratios^]. Each coefficient of this series can be calculated 



analytically by means of the asymptotic expansion algorithm described in [12 



To keep control of gauge invariance we work in a gauge with three different gauge 



parameters t^w, £z an d £ 7 . The corresponding free vector boson propagators Q2.4Q thus 



8 DIANA generates additional information, e.g. identifying symbols for the particles of the diagram and 
their masses, distribution of integration momenta, number of fermion loops etc. 

9 For diagrams with several different masses, there may exist several small parameters. In this case we 
apply different asymptotic expansions (see |28| ) one after another. 
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n< 1 > - O + O + Th 




exhibit the new masses \f^w m w an d \JXz m z in the propagators. This complicates the 
calculation enormously both in Tarasov's algorithm as well as in the asymptotic expansion 
approach. With the first method the presence of new masses in both the reduction formulae 
and the Gram determinants leads to cumbersome expressions which are difficult to simplify. 
In the asymptotic expansion approach the unphysical parameters £,w m w an< ^ £z m z define 
two new scales. Of course, in order to keep things manageable, we have to keep the number 
of scales as small as possible. This can be done by expanding diagrams about some fixed 
values of the gauge parameters. Three different regimes of expansion are feasible: a) £ — > oo, 
b) £ — > and c) £ — * 1. We choose the last possibility expanding the original propagators 
at £j = 1 in a Taylor series. For the purpose of checking the gauge invariance of our results 
it is sufficient to keep the first three terms of the expansion, so that the propagators of the 
vector bosons and associated Higgs scalar ghosts look like 

DUP) = -r^r [-9r + (1 - tv)^* - (! ~ fr) 2 ( f VP % 2 +■■), 
^ p 2 — my \ p 2 — rriy (p 2 — m v ) 2 J 

Av(p) = ^^fl-(l-^)V^ + ( 1 -^) 2 fV^V + ---V ( 3 - 8 ) 

P 2 - Tfly \ p 2 - my \p 2 - my J I 

where V = W, Z and the dots stand for the terms of higher order in (1 — £y), which we 
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don't take into account. At the same time we do not have any problems with the photon 
propagator and use it in its usual form. 

One drawback of our choice of the expansion about £j = 1 is that possible unphysical 
thresholds do not show up in our expansion and hence will not produce the correct imaginary 
part for a given diagram. Examples are the thresholds p 2 = 4£zM% (Z — > <p°<j) production, 
<ft° the neutral Higgs ghost) which is below the Z mass-shell p 2 = M| when £z < | or 
p 2 = ^y/M'y V (W ± — > ± 7 production, ^ the charged Higgs ghosts) which is below the W 
mass-shell p 2 = when £w < 1- in such cases, the above expansion does not preserve the 
analytical properties of diagrams exhibiting ghost particles. However, in any case unphysical 
degrees of freedom should not contribute to the physical width. While individual diagrams 
exhibit an imaginary part, in the sum of all diagrams it has to cancel. This cancellation 
is a consequence of the Slavnov- Taylor identities, which tell us how Higgs ghost, Faddeev- 
Popov ghosts and scalar components contained in the gauge boson fields decouple from 
observables like the physical width. Therefore, the above expansion indeed reproduces the 
correct gauge invariant result. At the one- loop level one can do an exact analytic calculation 
for the gauge with an arbitrary value of the gauge parameter. This has been done long 
time ago [HJ and, not surprisingly, the on-shell self-energies of the W- and Z-bosons are 
gauge invariant and do not exhibit any unphysical threshold (possible problems related to 
unphysical thresholds in the non-gauge invariant wave-function renormalization factor are 
also discussed in Hl6|). For the two-loop contribution, as we should, we get a gauge invariant 
on-shell limit which is real. Gauge cancellations are highly non-trivial and only happen if 
one is doing a perfectly consistent calculation. We have verified that for gauge parameters 
£ > 1 the imaginary part of the W and Z self-energies in the bosonic sector up to two- 
loops is zero, by applying the Cutkowsky rules and inspecting all possible two and three 
particle intermediate states allowed by the SM Lagrangian. While for £ > 1 the imaginary 
part is zero for each individual diagram, for small enough values of the gauge parameters a 
non-trivial cancellation must take place. An independent direct check of this is possible by 
considering the problem in the limit £ — > 0, for example. 

After we have made the expansion ( p.8|) , the bosonic contributions we are considering 
only depend on the three different masses m' 2 v ,m 2 z and m 2 H . One natural small parameter 
is the weak mixing parameter sin 2 9w = 1 — iriw/ m z < 0.25. We expand in this parameter 
and get rid in this way of (or mz)- For diagrams which contain Higgs boson linesQ we 
apply an asymptotic expansion with respect to a heavy Higgs mass. Taking into account 
the most recent lower bound on the value of the Higgs boson mass we are dealing with an 
expansion parameter m|/m 2 ^<0.64. This implies that we have to calculate quite a number of 
coefficient in the expansion in order to get a convincing result. We should mention that these 
expansions are well behaved because p 2 = M v in any case is below the possible thresholds. 
The convergence of such expansions can be easily checked for the one-loop case where the 
exact analytic result is available. 

In our case we find four different prototype structures with Higgs lines inside of the loops. 
The large mass expansion has been performed with the help of the packages TLAMM 
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We have 280 and 357 two-loop one-particle irreducible diagrams for the Z- and VF-boson, respectively. 
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and [51J once for Euclidean [3T| and once for Minkowski space-time, respectively. The 
corresponding topologies and set of subgraphs are given in Fig. |2|. The diagrams without 




Figure 2: The prototype diagrams and their subgraphs contributing to the large mass 
expansion for two-loop diagrams with heavy propagators. Thick and thin lines correspond 
to the heavy- and light-mass (massless) particle propagators, respectively. Dotted lines 
indicate the lines omitted in the subgraph. 



Higgs^ are nothing but single scale massive diagrams, when all internal masses are equal 
to the external momentum or zero. Such diagrams can be calculated analytically. For this 
purpose we use the packages ONSHELL2 j32| and another one written by O.V. |33| for the 
calculation of the set of the master integrals given in We find that only the following 
four prototypes are required (in terms of the notation used in [01): Flllll, F11110, 
F01111 and F01101. 

For independent verification of the input, the Feynman rules and the evaluation we 
performed calculations independently in Euclidean (M.Yu.K) and Minkowski (O.V.) space- 
time and got full agreement between them 0. 



11 At the two-loop level we find 336 one-particle irreducible diagrams without a Higgs line for the Z-boson 
and 435 for the T^-boson. 

12 The packages used in the c alculations and the results can be found at the following URL address 
ittp : //www-zeuthen . desy . de/ ~kalmykov/pole/pole .html. 
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4 UV renormalization in the MS scheme 



Here we describe in more detail the renormalization procedure. It is well known that the 
pole mass in QED/QCD is a gauge independent and infrared stable quantity to all orders 
of the loop expansion To renormalize the pole mass at the two-loop level requires 

to calculate the one-loop renormalization constants for all physical parameters (charge and 
masses), and the two- loop renormalization constant only for the mass itself. Not needed 
are the wave-function renormalizations or ghost (unphysical) sector renormalizations. The 
above mentioned basic properties of the pole mass are valid also in the SM. In order to obtain 
a gauge invariant result in the SM, however, we have to add in a proper way the tadpole 



contributions Hlq] . The tadpole terms are due to the vacuum expectation value (VEV) of 
the Higgs field, which does not vanish automatically. By a constant shift we can adjust the 
Higgs field to have vanishing VEV, however. Since here the Higgs field is integrated out in 
the path integral the result cannot depend on whether we perform a shift or not. Indeed, 
after we take into account all diagrams shown in Fig. [I] we find a gauge invariant result for 
the pole mass up to two loops in terms of the bare parameters. The tadpole contribution 
can be calculated either from tadpole diagrams^ or from Ward identities, which connect the 
tadpoles with the one-particle irreducible self-energies of the pseudo-Goldstone bosons 
at zero momentum 

nS I (0) + T = 0, (4.1) 

where <fi = + , <p~ , 0°. We performed both types of calculations and obtained full agreement. 
Diagrams contributing to the Z-boson pole mass do not contain infrared singularities at the 
two- loop level. Infrared fmiteness of the VT-boson mass was proven in |36[ . We also give an 
alternative proof of this statement for our case. 



In our calculation dimensional regularization |37] is used, which allows one to regularize 
both UV and IR singularities by the same parameter e related to the dimension of space-time 
by d = 4 — 2e. We first perform the UV-renormalization within the MS scheme in order to 
obtain finite results. In a next step we find the relation between the pole- and MS parameters. 
We adopt the convention that the MS-parameters are defined by multiplying each L-loop 
integral by the factor (exp(7)/47r) £L . Each loop picks an additional factor 1/167T 2 . 

4.1 One- loop charge renormalization 

The bare charge e and the MS charge^ e are related via eo = fi e e (l + Z^/e + 0(e 4 )), 

with the appropriate constant Z^. At the same time, the relation between the MS and the 
on-shell charge reads 

(i) 



e = e os (1 + 4s + °( e os 

where eos is defined by the Thompson limit of Compton scattering. The electromagnetic 
Ward-Takahashi identity implies that some of the diagrams cancel, such that Zq\ at the 
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These are the two-loop bubble diagrams for the process H in terms of QGRAF notation. 

MS 



14 All MS-parameters, like e, g, g' and all renormalization constants, like Z^s or zqs are /i-dependent 



quantities. 
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one-loop level can be written in terms of self-energies only |T5 



UV-fmiteness of ^Qg implies 



z w - -InW'fo) + sin6w n ^ (0) - - Z {1) (A 2) 

Zos ~ 2 ii ^ W + cos^ Mf e***' W 

(i) 



<°=H i -\\4? )■ (4 - 3) 



This may be confirmed also by using a renormalization group analysis of the SM keeping 
the Yang-Mills and the Higgs sector only. From the relation 



"2 = -2 + ~f2 ' ( 4 - 4 ) 
e <? <7 

where g' and g are the U(l) and the SU(2) gauge coupling constants, respectively, it is easy 
to deduce that 

The /3-functions may be calculated in the unbroken theory. They have been calculated 
in |38[| . We see that the above result is in agreement with (|4.3| ) if we take into account that 
(fj^d/dfi 2 ) e = — |e + j3 e . 



4.2 Mass renormalization 

We introduce the following notation for the mass renormalization constants 

mVo - <w ( t + M 4 m> + _^o_^> + ^4->) , ( 4,) 

where V stands for any of the bosons Z, W or if. In addition to the masses we have 
one coupling constant as a free parameter of the SM which we have chosen above to be 



e = gsinOyy. The one-loop mass counter-terms are well known [ I6[ . For the purely bosonic 
contributions we have 



(i i) 3 3 m| 3 m 2 H 

Zh = "2~4<: + 4<;' (4 - 7) 



_ 3/7% _ m 2 w _ 3 m| 3 m\ _ 17 

- 7~2 _2 2 _2 + /i _2 o' V^ 



4tt% r?% 2m l H m l w 4tt% 3 



4 lll} = -2^-3^-2-^ + ^4.-7^ + 1. (4-9) 
4 77% 77% 2 77%r?% 12 77% 6 



All masses here are MS masses and depend on the renormalization scale [i: rriy = m v (fi). It 
should be noted that unlike in the case of the couplings the mass renormalization constants 
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cannot be calculated from the unbroken gauge theory. Here in any case the calculation of 
the Feynman diagrams in the Standard Model is required. 

Let us comment about the somewhat unusual looking dependence of our MS renormal- 
ization constants on the particle masses. Since masses are induced by the Higgs mechanism 
we have the mass coupling relations 

m w =9^, rn z = \j g 2 + g' 2 | , m H = V2X v 

where v is the Higgs VEV and A the Higgs self-coupling. One peculiarity of the "spontaneous 
symmetry breaking" and the related mixing of states, like the 7 — Z mixing, leads to the 
non-polynomial nature of the perturbation expansion in the SM. Due to the mixing, the 



actual Z coupling reads y g 2 + g' 2 = gj cos^ etc. In the dimensionless mass ratios the 
factors v 2 drop out and we have in fact just ratios of couplings. To a large extent this is a 
trivial consequence of factorizing out powers of g 2 which cancels against such factors which 
appear in the denominators of the Zy^s. However, there are also true inverse powers of 
the Higgs self-coupling present. They originate in the tadpoles which we need taking into 
account for the sake of gauge invariance. 

Our results for the two-loop mass renormalization constants are as follows: 



£■(2,1) 63 mjj 3 m 2 H ?>m 2 H m 2 z 301 m\ 17 m| 53 
w ~ 64m^~4m^^8 m 4 w ~ 192^4 + 12m^ ~ T 



31 mi 17 ml 176 ml 59 m% 
_| £ ± YL _| £ — ^ (4.10) 

12 m^m^ 2 m 2 H 3 m 2 H 24m 2 H m\ v ) 
f2 . 2) 9 m% 13 nr n 65 m% 35 my 334 27 m% 

Zw — 7777 a I 77 n \~ 7777 a 77 n I 77 I 



' w 32m$ v 8 m\r 32mfy 8 m\r 9 4 m\m\[ 

„ ml n t ml 5 m% „ mi? „ mi 9 m% , , , , s 

+6^ + 34-^ - -^r-V + 9 ^ + 9^ + 7^^, 4.11 
m # m # ^ m H m w m H m H ^ mjjm^y 



£.(2,1) 63 mfj 3 m^ 3m 2 H m 2 z 11 7 m| 64 m^ 253 mf 
z 64 m^ 4 m^ 8 m^ 3 6 m^ 3 m| 192 m^ 

176 m^ 17 m| 31 m| | 59 m| 
3 m^ 2 m^ 12m 2 H m 2 v 2km\m\ v ' 1 
(2,2) 9 m^ lm 2 H lm 2 H m 2 z 21m 2 H 55 11 m| 629 m| 

z ~ ~32 ^4 + 4^ ~ 8 m^ + ~ IF + 12^ + 288 m^ 



11 m| 



ozm w 4 m w a m w t m z u 

245 m^ 27m^ m| m^ 7 m| 

6 2 ' r> 1 ' 1 ' 11 nil A 1 A 

m z 2 m H m H m z m H 2 m^m^ 4 m H m^ 

m%j mi 9 ml 
+9 _»l + 9_^ + 4.13 

Again the higher order pole terms 1/e 2 can be checked by means of the appropriate 

renormalization group equation. Let us write the relation m vo = my (l + J2 n Zy^/e 71 ^ 
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which connects bare and renormalized masses and introduces the anomalous dimension of 
the mass 7y = (/i 2 d/d^i 2 ) mm 2 /. Taking into account that (n 2 d/dn 2 )my = and repeating 
the calculations given in |39| we find 

TV = i 9 |4", (4.14) 

For each function like 7y or Zy^ we may perform a loop expansion 7y = Z)^ =1 (5 ,2 /167r 2 ) fc 7y ^ 
or Zy^ = J2fc=n(9 2 /ISvr 2 )^ Zy' kS) and ( f4.15[ ) can be written for each loop correction separately. 
In particular, for n = 1, we have 

7?> = 4 W) , 

such that the coefficient of the n = 2 poles can be checked via 

(4">) 2 + + E *8Mjatf* = 24 2 ' • (4.16) 

y y i L " n i 

The value of = — f§y§^2 ma y be calculated from the relation 

om • Q 1 cos 2 # w ( 9 2 \ ( (l) . fl m 



(4.17) 



where (3^ is given in ( |4.5| ), and cos 2 9w = ^w/rnz- 

An independent verification of the 1/e terms can be obtained from the relationship 
e 2 = g 2 sin 2 9\y, which is valid for bare and MS renormalized quantities. Differentiating 
the renormalized quantities with respect to In /i 2 we find 

e 4 + &\ = gft g sin 2 9 W - \g 2 { lw - lz ) cos 2 6 W , (4.18) 

or for the two-loop case 

g 2 ^ sin 4 dw-^f- sin 2 W cos 2 9 W = - (J^J ' (z^ - zf x) ) cos 2 9 W . (4.19) 



The two-loop /3-functions for g and g' are given in |40| and read 

1 g' 3 , 1 g' 5 , 3 g' 3 g 2 



ft 



" ~ 1216tt 2 4(16tt 2 ) 2 4(16vr 2 ) 2 ' 
43 g 3 259 g 5 1 g 3 g' 2 
UW^ 2 ~ 12 (16tt 2 ) 2 + 4 (167T 2 ) 1 



= -^J* *u*_y x y y , v 



12 



The fact that after UV renormalization we get a finite result confirms the infrared finite- 
ness of the bosonic contribution to the pole mass. 

Let us write now the RG equation for the effective Fermi constant Gp. Gf is usually 
defined as a low energy constant in one-to-one correspondence with the muon lifetime. How- 
ever, if we consider physics at higher energies a parameterization in terms of low energy 
constants may lead to large radiative corrections. Much in the same way as the fine struc- 
ture constant a often is replaced by the effective running fine structure constant a(fi) we 
expect that Gf should be replaced by an effective version of it at higher energies. Unlike in 
the case of a, however, because of the smallness of the light fermion Yukawa couplings, Gf 
starts to run effectively only at scales beyond the VT-pair production threshold (see below). 
The bare relations are 

<o = ^ (4-21) 



\ = V2G F , . (4.22) 



Both these relations are valid also for MS renormalized parameters, so that their differenti- 
ation with respect to In /i 2 gives rise to the relation 

1w = 2^- 1Gf , (4.23) 
9 

where we introduced the anomalous dimension of the Fermi constant 7g f = (/x 2 dj dfj, 2 ) In Gf- 
Its loop expansion looks like 7c F = 7g f + lei + ' ' '• At the one- loop level we have 

7g» = 




i(2m i w + m i z )+m i H -AY, 



m 



f 



m 



-2£mJ (4.24) 



and we used the results of Ref. for the fermion contributions. The first term proportional 
1 / m 2 H is the contribution from the tadpoles. The appearance of the tadpole terms is some- 
what mysterious, since we know that in renormalized observables tadpoles drop out. Here 
they seem to contribute to the renormalization group evolution of the Fermi constant. In 
any case the tadpoles are present in the relationship between the bare and the renormalized 
parameters. At the two-loop level, our results allows us to write the bosonic corrections only. 
They are given by 

(2), boson 

1g f — A ~ AZj w 
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2g 4 \ 63 m\j 3 m\ 3 m\m\ 301 m% 7 m\ 25 

i 77 t r 



'167T 2 ) 2 [ 64 myy 4:771^ 8 192 677^ ' 6 

,31 m| _ 17 m| _ 176^ + 59 m| j 
12 m^m 2 ^ 2 m 2 ^ 3 m 2 ^ 24 rn? H m'\ v J 



The equations fl42jj ) and ( ggg ) are written in MS scheme. As usual in this scheme, in 



solving the renormalization group equation 

G F (fi) = G F (fi) 1g f 

the decoupling of the heavy particles has to be performed "by hand". This means that 
for low values of the energy scale /i, when \i < m H ,m w ,m z , the bosonic terms on the 
r.h.s. are equal to zero while the light fermion contributions proportional to G F rn?^ are tiny. 
Consequently, below the W mass, the effective Fermi constant does practically not change 
with scale. Obviously, the running of Gf only starts at about ii ~ m?, when the scale of a 
process exceeds the masses of the bosons. Also the top quark will contribute once we have 
passed its threshold. 

5 Results, discussion and conclusion 

After UV renormalization the pole mass (see ( |2.3j )) 

My = ml + n« + up + n«ng>' (5.1) 

is represented in terms of finite MS renormalized amplitudes 



n« = IT?V,m 2 



Vi 



, F.P. 



F.P. denotes the MS finite part prescription. The calculation of the one-loop MS renormal- 
ized amplitude is well known. We get it by rewriting the bare expression in terms of MS 
parameters according to ( |4.6|) and ( |C.1|) 

ftj? = ^ [rnl v -m v + ml v ^X^) 

= rn v (fi) lim (-Z^ + ig) = m^) 2 &% a X® ■ (5.2) 

167r^ sm W e^o \£ ' / lxyK z surVw 

We restrict ourselves to a consideration of the bosonic part X^y oson . The renormaliza- 
tion of the off-shell self-energy functions IIq 1 ^' and Il^y is more complicated. Besides the 
renormalization of the physical parameters, it would require us to perform order-by-order, 
the wave-function renormalization as well as the renormalization of the ghost sector, in 
particular of the unphysical gauge parameters. However, the MS renormalization of the 
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combination H^ v + H^yH^y' it much simpler. The subtraction of sub-divergencies in this 
case is reduced to the one-loop renormalization of the charge and the physical masses only, 
while the wave-function renormalization or the renormalization of the unphysical sector is 
not needed. Accordingly, as a genuine two-loop counter-term, only the mass renormalization 
occurs. 

The full two-loop MS renormalized amplitude can be written in the form 



r(!) tt(!) 



fr(2) , frWfrW' = J n (2) + n (1) TT (1) ' 

lly -t- lly lly _ < ll y + IIq yllgy 



+m v (n) 



1 



+m v (n) 



2 r 



£ V 16tt 2 sin 2 6 



w , 



+ 



MS 

As 2 



lim [ rr (2) 4-n (1) n (1)/ 

e— >0 y 



2 



<9 



(9m 2 



167r 2 sin 2 6 



w 



7 (2,1) , 1 7 (2,2) 



(i) 



(5.3) 



where the sum runs over all species of particles j 



Z, W, H and 



As 2 



cos 2 # 



sin 2 6 W 



7(1,1) 



7 sin 2 . 



The functions Z^X^ and are defined in Q-(p|), (p^)-(gl3|), (Q and QCp 



-(I) 



respectively. Again, throughout, we take into account only the bosonic corrections. The 
r.h.s. of ( |5.3|) is given by the bare two-loop contributions, the two- loop contribution ob- 
tained by expanding the bare parameters in the one-loop amplitude (subtraction of the 
sub-divergences) and the genuine two-loop subtractions. The second type of contribution 



follows from writing IT 1 ) 



m; 



(it?) ^0 



(1) 



and by utilizing 

(i) 



~M. We would like to stress, 



"0 V167T 

that each of the one-loop bare amplitudes and Hq 1 ^' has to be expanded up to linear 
terms in e. Together with the singular \ je terms the contributions linear in e yield additional 
finite contributions in the limit e — * 0. In the product n^yllQ 1 ]/' the parameters may be 
identified by the MS renormalized quantities, but all functions, Aq and Bo have to be taken 
in d = 4 — 2e dimensions and must be expanded up to terms linear in e (see Appendix A). 

As mentioned earlier, for the purely bosonic contributions alone the imaginary part of 
Il(p 2 ) on the mass-shell is zero at the two-loop level. This is due to the fact that in the 
bosonic sector we have the physical masses m 7 = 0, Mz, My/ and M# and by inspection of 
the possible two and three particle intermediate states one observes that all physical thresh- 
olds lie above the mass shells of the W and Z bosons, i.e., the self-energies of the massive 
gauge bosons develop an imaginary part only at p 2 > My (to two loops in the SM). On kine- 
matical grounds imaginary parts could show up from the Higgs or Faddeev-Popov ghosts, 
which have square masses £yM v , for small values of the gauge parameter. However, as we 
have verified, the two- loop on-shell self-energies are gauge independent. This implies that 
ghost contributions have to cancel and hence cannot contribute to the imaginary part. Thus 
sp = My in our case. In higher orders for the Z-propagator one gets an imaginary part as 
soon as p 2 > 0, from diagrams like 
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7 




7 



For the VT-propagator an imaginary part is only possible for p 2 > M^, because charge 
conservation requires at least one W in any physical intermediate state. 
For the massive gauge bosons Z and W we write 



My 



1 + 



rn 



v 



lGir 2 sin 6 W 



167r 2 sin 2 6 



X 



(2) 



w 



(5.4) 



where both e and sinOw are to be taken in the MS scheme. 

The one-loop coefficients Xy^ for Z, W and H are known of course as exact results. 

(2) 

We write them down for completeness in Appendix B. For the coefficients X v we make an 
expansion and perform them as double series: in sin 2 9yy and in the mass ratio m v /m 2 H . 
We have calculated the first six terms of the expansion with respect to sin 2 6w and the first 
six terms with respect to mass ratio m v /m 2 H . The analytical values of these coefficients are 
presented in Appendix D. These represent our main result. 

Sometimes in massive multi-loop calculations the so-called modified MS scheme (MMS) 
is used The difference between MS and MMS is that in the former scheme each loop 
is multiplied by (e 7 /47r) e while in the latter the normalization factor is l/(47r) e /r(l + e), 
which yields a difference at the two-loop level. It has been shown that in QCD both schemes 
reproduce the same formula for the mass relation analogous to (|5.1| ) ||42|| . We have checked 
that the same holds true for the pole masses in the Standard Model. 

Very often the inverse of ( |5.1| ) is required. To that end we have to express all MS 
parameters in terms of on-shell ones. Thus 



v 



n 



(i) 



n (2) , n (1) n (1)/ 

LL o,v + ii o,v ii o,y 



MS 



-B^><"^>-(Ae)<'>^> 



m 3 = M j > e=e OS 



(5.5) 



where the sum runs over all species of particles j = Z, W, H and 



(Am 2 )« = -Refi^ 



-Ml 



'OS 



mj=Mj, e=eos 



IQtt 2 sin 6w 



X 



(1) 

V 



m 2 =M 2 
j J 



stands for the self-energy of the jth particle at p 2 = m 2 in the MS scheme and parameters 
replaced by the on-shell ones. Note that in the above relation we had to perform a change 
from the MS to the on-shell scheme also for the electric charge 



e{n 2 ) = e os 




os ,_ ln 



(5.6) 
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with eQ S /47r = a ~ 1/137. Accordingly, since IT 1 ) depends on e by an overall factor e 2 only, 



7 In 



. ^ 2 . 



n 



(i) 



By identifying m 2 



v 



m 



wo 



My + SM V ( |5.5|) is the relationship appropriate to obtain 



the on-shell gauge-boson mass counter-terms SM V : 



5My 



-Re 



tt(i) , fr(2) , frWrrW' 
il y,o + ^v.o + ■ u v;o li v,o 



2^(1) 



9 



3 

■ ^OS 



dm 2 ny '° 



'j',0 



1 M 2 



<9 



(i) 



(fc)«"-n, 



m ],0= M j >eo=eos 



(5.7) 



in terms of the original bare on-shell amplitudes 



fr(<) 



y,05 



p 2 =M v ,m? i0 =M?, e =e O s 



and the bare on-shell counter-terms <5M 2 and <5e. The second equality gives <5My in terms 
of the singular factor Zy^ = m vo /m v (fi) given in ( |4.6| ) and the finite factor Zqs = m v /M v 
given in (|5.5|) . These will be needed in two- loop calculations of observables in the on-shell 
scheme. 

We now turn to a discussion of our results which we obtained for the relationship fl5.5| ). 
For our numerical analysis we used the following values of the pole masses: M\y = 80.419 
GeV and Mz = 91.188 GeV and a = 1/137.036. We first investigate numerically the re- 
lationship between the MS-mass my(ju) and the pole-mass My. In Figs. |] and | we plot 
A z = m\{M z ) /Ml - 1 as function of the Hig gs mass Mh for intermediate and heavier Higgs 
masses, respectively. For the one-loop corrections the exact analytical functions are evalu- 
ated while for the two-loop results we utilize all coefficients of our expansion. Analogously, 
in Figs. || and | the Higgs mass dependence of = m 2 y (M M /)/M^ / — 1 is depicted for the 
same ranges of the Higgs mass. As we can see, for a "light" Higgs of mass less than about 
200 GeV the two-loop corrections are small as compared to the one-loop ones. However, at 
a Higgs mass of about 220 GeV the absolute value of the two-loop correction is of the same 
size as the one-loop result, such that the two-loop corrections start to play an essential role. 

Our analysis shows (see Figs. £| and ^) that the perturbative expansion looses its meaning 
for large Higgs masses (strong coupling regime). In the relationship between the MS - and 
the pole-masses corrections exceed the 50% level around 880 GeV (for /i = My). The size 
of the corrections depend on the choice of the renormalization scale [i. 

Since our results have been obtained by an expansion in m v /m 2 H (i.e., for a heavy Higgs), 
one of the main questions which remains to be considered is the validity of our results as 
we approach lighter Higgs masses. Note that we are dealing with an asymptotic expansion 
(which means the radius of convergence is zero) and thus we can answer the question about 
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where the expansion is reliable only empirically; at least as long as an exact result is not 
available for a direct comparison. Since the structure of the one- and two-loop corrections 
for both massive gauge bosons is very similar, we are going to investigate in the following 
the problems of convergence of our results for the Z-boson only. For the IV-boson the 
convergence is better because of the smaller value of the V^-mass. Obviously our series- 
expansion breaks down for a "light" Higgs, when m#>m^. Unexpectedly, we find that the 
two-loop corrections remain numerically small for values of the Higgs mass even down to 
100 GeV (below about 0.2%). However, we observe a steep raise of the correction which 
signals that the expansion becomes unreliable below about 130 GeV0. We also should keep 
in mind that, on the level of the size of the two- loop corrections, which are very small in the 
region around 150 GeV for fi ~ m z , the latter depend substantially on the choice of the MS 
renormalization-scale /i, which causes an essentially constant shift in the quantities shown 
in the Figures which follow. 

In Fig. [7| we show the dependences of the two-loop corrections to A z = m z {M z )/Ml — 
1 as a function of the Higgs mass and the number of coefficients of the expansion used 
for their evaluation. For a "light" Higgs the difference between the result, including six 
coefficients of the expansion, and the results obtained by including only the first three of 
them (leading, next-to-leading and next-to-next-to-leading) are numerically small. However, 
the convergence of the series is not satisfactory. The coefficients grow fast beyond the first 
four terms and the correction starts growing fast. For a heavy Higgs with a mass of more than 
300 GeV the convergence is much better, so that we omit the corresponding plot. We only 
mention, that there is no essential differences between the full result and the next-to-leading 
one. Similarly, the dependence of the two-loop corrections on the number of coefficients of 
the expansion with respect to sin 2 9w for a "light" Higgs is illustrated in Fig. H Here unlike 
in the previous "light" case, the first coefficient is relatively large already, while the higher 
terms of the expansion do not alter the result in a significant manner. Again, we observe 
that results become unreliable below about 130 GeV. 

Finally we analyze the Higgs mass dependence of sin 2 6w The relation between the MS 
weak mixing parameter and its version in terms of the pole masses reads 

sm9w ~ mi~ Ta [t^w^w 



i 



L z 

Ml 

Ml) Ml 



where we adopted the notation, m v /My = 1 + 8y + 8y . It turns out that in (|5\8| ) all Mfj 
terms cancel. The corresponding results are depicted in Figs. and [TO. Again for a "light" 



Higgs the two-loop corrections are small, while for a heavier Higgs particle the corrections 
become large. Again, below about 120 GeV we cannot trust the series-expansion any longer. 



15 We have checked the convergence of the corresponding expansion at the one-loop level, where the exact 
result is available. Deviations start to show up below about 145 GeV and grow to about 30% at 100 GeV. 
The best approximation in this case is based on the first 3 or 4 coefficients. In this case, beyond the first 
few terms, the series starts to diverge in a way typical for an asymptotic expansion 
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We would like to mention that the presence of mjj corrections in the relation 
does not contradict Veltman's screening theorem [43 which states, that the L-loop Higgs 
dependence of a physical observable is at most of the form {jn 2 H ) L ~ l ln L m 2 H for large Higgs 
masses. This theorem applies to physical observables like cross sections and asymmetries, 
whereas our formula is nothing but a relation between parameters of two different schemes. 

In conclusion, the main results of our paper are the following: (i) we have presented 
an independent proof of gauge invariance and infrared stability of the 2-loop electroweak 
bosonic corrections to the pole of the gauge boson propagators; (ii) analytical results for a 
number of coefficients of the expansion in sin 2 9w = 1 — m 'wl m 'z and my/m 2 H for the 2-loop 
electroweak bosonic corrections are given for the on-shell mass counter-terms ( |5.7| ) and the 
relationship between MS and on-shell masses of the gauge bosons W and Z. All calculations 
have been performed in the electroweak Standard Model. 



Note added: After completion of our paper we received the preprint ||52[| , which presents 
a complete SM calculation of the Higgs mass dependent terms to the observable Ar which 
determines the Mw — Mz relationship given a, and Mz as input parameters. 
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A The one-loop master integral and its e-expansion 

For the two-loop calculation we have to take into account the part proportional to e of the 
one-loop propagator type integral 16 

J= f ^ (A1) 

J {q 2 - ml + i0){{k - q) 2 - m 2 2 + iO) ' K ' } 

where d = 4 — 2e and "+i0" is the causal prescription for the propagator. Its finite part has 
been presented in ||44j| , the 0(e)-part in | 4^| , the terms up to order e 3 can be extracted by 



means of Eq. (A. 3) of |H| and an all order e-expansion was obtained in ||46|| . We write the 
part of J linear in e in a form suitable for the implementation in FORM 



J 



17T 



2(1 -2e) 



m. 



-2e 



-m, 



-2c 



+ 



m\ 



-mt. 



e 



e k 2 



m 



-2s 



—m. 



-2e 



16 In II it is denoted as J (2) (4 - 2e; 1, 1). 

17 The higher order e terms can be extracted from 
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k 2 



arccos 



m 2 + m\ — k T 
2m\m 2 



1-eln 



-A(m 2 , mf, /c 2 ) 
fc 2 



+2 £ (C1 2 (n) - Cl 2 (tt - n) + Cl 2 (r 2 ) - Cl 2 (tt - r 2 ; 



+ 0(e 2 



(A.2) 



where X^m^m 2 ,, k 2 ) = [m\ + m\ + fc 4 — 2m\k 2 — 2m 2 ,k 2 — 2mfm 2 ,) and the angles T{ are 
defined (see via 



COS T\ 



k 2 — m\ + m 2 



2mo\/k 2 



cosr 2 



fc 2 + m 2 — m 2 



2m\\/k 2 



(A.3) 



C1 2 (0) is the Clausen function Cl 2 (0) = | [Li 2 (e ie ) -Li 2 ( e - ie )]. The expansion Q is 
directly applicable in the region where A < 0, i.e. when (mi — m 2 ) 2 < k 2 < (mi +m 2 ) 2 . For 



the region A > we need the proper analytic continuation which has been given in ||47|| . Let 
us briefly describe it here. First of all, we rewrite ( |A.2| ) in the form (see section (2.2) of |47j 
for details) 



2 _ £ r(l+e) /mr 2£ +m 2 



J = m 



-2e 



2(1 -2e) 



m\—m\ 



m 1 2e —m 2 



2s 



-\(m 2 ,m 2 ,k 2 )} 1/2 ~ £ f 2r 2 (l - e) 



— l- 



(i - pi - p 2 ) + - E I p*(-^r e - (i - p*)(-^) £ 



+2£Efp»(-^)" Li 



\ e r(l-2e) 
2 { Zi ) - (1 - Pi )(-^) £ Li 2 (lM) + O(e) 



(A.4) 



where pi and p 2 are some numbers, which we will define later and 

2 



\{m\ , m 2 , fc 2 ) + m 2 — m\ — /c 5 



4m 2 /c 2 



2 2 



A(mf , m 2 , A; 2 ) — m 2 + m\ — k" 



Am 2 k 2 

Firstly, we note that the causal prescription amounts to the following rule for A (A > 0) 

ln(— A(m 2 , m 2 , k 2 )) = ln(A(m 2 , m 2 ., k 2 )) — iir, 



(A.5) 



\/-A(mf,mi,fc 2 ) = -iy / A(mf,m|, fc 2 ). 



The function Li 2 (z) is real for real z and |z| < 1. For real z and \z\ > 1 we change argument 
z — > I/2 using the relation Q 



Li 2 (z)+Li 2 Q =-^ln 2 (-^)-C 2 
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For the higher order poly-logarithm the relation is 



Li„(z) + (-l)"Li„ ( - ) = -i ln"(-z) - £ 



/ 2 1 1 n—1r 1 \ 
In (-z) 



^ (n-2r)! 



<2 
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by which an imaginary part shows up. This change of variables can be done from the very 
beginning in (|A.4j) by an appropriate choice of the values of the coefficients pf 



< Zj < 1 =>■ pj — 1; ln(— Zj) = hi(zj) + m, 
Zj > 1 =>• pj = 0; hx(—Zj) = hi(zj) — yk. 

Assuming m\ < m 2 in the following we have 

• for k 2 < (mi — m 2 ) 2 =>- Z\ < 1, z 2 > 1 



. , F r(l + e) (m A 2e +m 2 2e m 2 -ml 
J = m 2 ' 



2(1 -2e) V £ 

J\{m\,ml, k 2 ) ( 
+ ^ p jhi(^ 2 ) + £ 

- ln(ziz 2 ) In 

for k 2 > (mi + m 2 ) 2 =^> z x < 1, z 2 < 1 



+ 



ek< 



m, 2s -m« 2e 



■2U, [ ^) - 2Li 2 (^i) - I In 2 zi + 1 In 2 z 2 



A(m 2 , m 2 , 



+ 0{e 2 ) 



(A.6) 



J = m 



r(l + e) /mr 2e +m 2 - 2e 



2(1 - 2e) I £ 



1 2 - f m x e - m 2 l£ 



ek 2 



J\{m\,ml, k 2 ) f 
+ - — S ln(ziZ2) + £ 



A- 2 



1 1 
-- In 2 zi - - In 2 z 2 - ln(ziz 2 ) In 



8C 2 -2Li 2 (zi) -2Li 2 (s 2 ) 
/ A(m 2 , m|, A; 2 )' 



fc 2 



+i7T 



2 - 2e In 



A(m 2 , m|, k 
k 2 



+ 0(e' 



(A.7) 



In a similar manner, starting from Eq. (2.17) of [17] and performing an analytical continua- 
tion, it is possible to obtain the higher order terms of the e expansion In particular, the 
imaginary part of J in each order of e coincides with that obtained from the exact result 



[0 



ImJ = m9 Ik — (mi + m 2 ) 



^A(m 2 ,m 2 ,P) /A(m 2 ,m 2 ,A; 2 )\- £ T(l - e) 



k 2 



k 2 



r(2-2e)' 



In the limit, when one of the masses vanishes, the result is [E7 



J \m 1= 0, m2 =m = ™ 2 £m 2£ 



r(l +e) [ 1 1-u 
(1 -2e) I e 2^T 



:i- M )- 2e -i 



:i-u 



,l-2e 



-eLi 2 ( M )+0(£ 2 )|, 
(A.8) 



in 



A collection of useful expressions for the one-loop two-point function is given also in Appendix A of [|49| 
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with u = k 2 /m 2 . 

The transition from the bare parameters to the renormalized ones requires differenti- 
ations of the one-loop propagators with respect to all parameters, couplings, masses and 
external momentum. The integrals obtained thereby can be reduces again to integrals of 



type ( |A.4j ) plus simpler bubble integrals. The expansion of the propagators with respect to 
small parameters (ratios of the masses or momenta and masses) can be extracted from the 



exact analytical results written in terms of hyper-geometric functions (see |50| ) 



B MS vs. pole masses at one-loop 



In this Appendix we present, for completeness, the well know [16|] one- loop relations between 
pole and MS masses of gauge bosons. Using the following notation 
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C Unrenormalized one-loop expressions in d dimension 

The computation of higher loop corrections requires a deeper expansion in e of lower order 
terms. In this Appendix we present for completeness the results for unrenormalized one-loop 
corrections to the relation between pole and MS masses of the gauge bosons in arbitrary 
dimension d without expanding it in e. Using the following notation 
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In above formulae we use the following functions 
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m Ui and m^ denote the masses of corresponding up- and down-quarks, N c is a number of 
color (N c = 3) and Kij is the element of the Kobayashi-Maskawa matrix. 



D MS vs. pole masses at two-loop 



After expansion of the diagrams with respect to sin 2 6y/ we get rid of one of the boson masses 
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and write the functions X v introduced in (|5.4| ) in the form 
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In particular for the Z boson propagator we eliminate mw and vice versa. Consequently, the 
coefficients AY in the above formula are functions of the Higgs mass and one of the boson 



25 



masses. We expand this function with respect to my/m 2 H 
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and calculate analytically the first six coefficients. This is not a naive Taylor expansion. The 
general rules for asymptotic expansions |T2j allow us to extract also logarithmic dependences, 
or in other words, to preserve all analytical properties of the original diagrams. In the result 
of the asymptotic expansion all propagator diagrams are reduced to single scale massive 
diagrams (including the two-loop bubbles). As a consequence, the finite as well as the e- 
part of the corresponding diagrams, are characterized by a restricted set of transcendental 
numbers |5T| which may appear in the coefficients Aij. We find the following constants: 
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Furthermore, ln(m 2 H ) denotes In (m' H / u 2 ) where \i is the 't Hooft scale. We also introduce 
the notation wh = rn^/mjj- and Zh = m 2 z /m 2 H . 



D.l Analytical results for the two-loop finite part of W-boson 



+ w H 



W-tSo + fS 2 + f|7r 2 -f H w H ) + § In(i^) ln(m|) 



519 
32 



Mm H ) + | Hm 2 H )S + f (Mm H ) 



+ w 



H 



8311573 153 c 243 c 87651 c 47893 _2 , 459 a/q\ 
20736 8 D 3 4 ^0 i 28 D 2 3456 ^+16 <,V>) 



399943 



1705 
16 

1595 
16 



\n{w H )S + ^ (Hw H ) f + iff ln(w H ) \n(m 2 H 



1728 

356585 i„f„.,2 
1728 



ln(wtf) 



\n(m H ) 



\n(m 2 H )S + 



10013 
288 



(\n(m 2 H ] 



+ w H 

517 



10730119 _ 3451 
86400 



3451 c 
-gg-^O 



:i " ,7 ' ,! ' ^ ^ ^ - ^H? Hwh) 
Hw H )S + {Hw H )Y + ^ ln(w H ) ln(m 2 H ] 



>2 
\2 



1080 
1743 



30383 
192 



ln(m^) 



+ 2 -f Hm 2 H )S + ^ (\n(m 2 H ] 



+ w 



H 



276774409 , 86473 c 3807 c i 859 _2 168691 i T1 /'„„ \ _j_ 10323 \ c 

1296000 + "3200" ^0 - ^2 + ^Tf VT - -j^ til{W H ) + -j^- til{W H )b 



26 



-W ( ln «) 2 + ^ 



\n(w H ) \n(m 2 H ) + ^ ln(m^) + f (ln(m^ 



13424129921 , 416057 o _ 452727 c i 1395311 2 _ 3016477 | / \ 
14112000 ' 2400 1120 2 ' 30240 10080 "H" 7 ^ 



7739 
40 



\n(w H )S + ^ Qn(w H )Y + ^ M^ir) ln(m£) + 



2 . 12347 



6333137 
30240 



ln(m|) 



i + i5 + ^7r 2 + | ln(^) - £ 1„(^) l„( m ^) _ | ln(m ^ _ 3j ln(m ^ 



w 



ii 



15491 I 65 c 11183 C i aa Q 4239 c i 8837 _2 22 2 i„/q\ ,67/ 



_ i4| ln( ^ } _ 209 Hwh)Sq - m (\n(w H )) 2 - f ln(^) ln(m|) + ff ln(m^) 



ln(m^ - £ (ln(m| 



4 



14221 , 977 o I 8181 o , 5705 _2 5923 \ , 77 l^A,,, \ c 



+ W (Hwh)) 2 + 1 ln(^) ln(m^) - f£ ln(m^) + f ln(m^)5 + 12 (ln(m^))' 



H 



43842631 , 259993 



i 259993 Q I 4329 c i 21149 ^2 , 26903 \ _i 4431 1ti /.„ \ o 



40579 
720 



(Hw H )) 2 + \n(w H ) \n(m 2 H ) + ff \n(m 2 H ) + 27 (ln(m^ 



367144853 , 2094431 c _i_ 

1 7onn «-?0 "T 



518400 



7200 



29151 
320 



So 



210541 2 | 8105093 



2880 



7T^ + 



43200 



ln(wtf) 



ln(^)5 + ™ (ln(^)) 2 - f ln(^) ln(m^) - ff ln(m 2 H ) 



■ 2 i + 3 i S ° + T^ 2 + 2 i Hwh) - Te Hwh) Hm 2 H ) + f Hm 2 H ) - f Hm 2 H )S 



w 



H 



439321 



+ ^S 3 + ^S -US 1 + 



14469 



So- 



19417 ^2 



3456 1 36 6 1 144 u i 64 ^ 576 

- if Cs - 1 Hwh) - m H™h)S - fg (ln(^)) 2 
+ f| ln(^) \n(m 2 H ) + ff ln«) - ? ln(m^)5 + § (ln(m^)) 



7r 2 -16 7r 2 ln(2) + ^7T 2 ln(3) 



?4 



304691 , 15601 



15601 Q i 10659 Q i 9167 J2 i 4267 \-.(.,, \ i 1073 l„^„ \C 



3456 



3509 
" 144 



(ln(^)) 2 + f In( WH ) ln«) + f In(m^) + x -f Hm 2 H )S + f (ln(ml)) 1 



27 



w% 



107428459 , 990193 c i 5103 c i 111 Jl , 10042841 w„„ \ i 10771 l / *\ c 

432000" + ~%0b~~ ^0 + ~W ^2 + — + 86400 m{w H ) + ln{w H )b 



+ 3||9 (l n (^))2 + 821 Hwh) Hw 2 h) + 1|1 ln(m 2^ + I ^ n(m 2^ 



H 



3726216829 , 2622653 c 1 98229 q 1 1160293 _2 , 21966593 l-A,,, \ 
"3600~~ ^° + ~320~ ^2 + 8640 ^ + 43200 ln \ W H) 



2592000 



+ 



27221 
60 



\2 , 1 



w H 



W 2 H 



-2 + f S + A ^ + I _ A i n (^) in ( m ^ + | ln(m ^ _ | ln(m 2^ o 

799 199 c 26975 c 5c, 4745 c 1 1669 _2 , 104 Jl \„(o\ 425 A/q\ 
-— - "54 ^3 - ^44- ^0 - 3 ^1 + — ^2 + — VT + ^7T ln(d) - -gg- C(JJ 



+ 2f ln(^) - 1h(«;h)5o - f (ln(^)) 2 + f lnfan) ln(m^) + f ln(rr&) 



7795 

" 288 



857575 1 196669 c , 49193 c 1 13715 Jl , 26065 w„„ \ , 2429 i „/,,,, No 
"69TF + 5l84- b + "W 6 2 + 7T + \n{w H ) + w lllfu^Po 

(ln(^)) 2 + ff \n{w H ) ln(m 2 H ) + ln(m 2 H ) + ±f ln(m^)5 + 8 (ln(m^)) ; 



3314383 , 1531019 
10800 



9000 



5-0 + 



42717 
160 



S 2 + 



53773 J2 , 6187879 



1440 



7T' + 



21600 



ln(^) + f£ Hw H )S 



fiflf (ln(^)) 2 + ^ ln(w H ) ln(m 2 H ) + ^ ln(m^) + 81 (\n(m 2 H 

- - ln(/r,/) 



3404318779 _j_ 1761697 c i 288607 c i 1098061 2 , 4020149 
1296 480 2 i 4320 71 "r 3600 



1296000 



38929 
54 



hi(w H )S t 



43325 

-r 



(ln( W „)) 2 + f ln(^) ln(m|) + ff ln(m 2 H ) 



3095 + if So + £ - 2 + I M«*) - & ^(^) ln«) + 1 ln«) 



1152 



144 



]n(m 2 H )S 



w 



H 



40045 841 q 747209 C i 1 C i 125399 c _i_ 660755 Jl , 68 Jl ]„fo\ 791 /■ 
'^32~ _ 216 ^3 - "518T" ^0 + 9 ^1 + -354- ^2 + ^HoT 7T + ^ 7T lll^j - m (,3 



-ff ln(^) - H In(^)5 - ff (ln(^)) 2 + f ln(^) ln(m|) + ff ln(m 2 H ) 



28 



H ln(m^ + f 



+ w 



H 



803113 _|_ 323699 
7776 



5184 



C , 33121 c I 28757 „2 , 85597 
"I 5fi~~ *- , 2 T 2592 ~ l ~ 



864 



ln(w H ) 



30247 
1296 



\n(w H )S 



+ H (Hw H )) 2 + f Ib(«, h ) In(m^) + ^ ln(rr&) + *f in(m |) <? + 2 ( ln(m | 



1919013 , 40790531 c i 
4000 o^oonn ^0 ~T 



869239 o i 1015361 ^2 , 16268087 w„„ \ , 7991 i„/.„ \o 
259200 -0 -T ^920" ^2 + 17280 71 + 28800 ln l W ^J + W "H^H Po 



+ ™ (ln(^)) 2 + + ^ + ¥ o«)y 



305951057 _i_ 414264923 
™ " 194400 



72000 



S + 



265819 
240 



S 2 + 



955459 Jl , 11862391 



2160 



7T^ + 



5400 



ln(w H ) 



3022481 
3240 



\n(w H )S + ^ (Hw H )f + ^ ln(^) ln(m^) + Sg* ln(m^) 



^5 



"US + W + ft vr 2 + § ln(^) - ft ln(^) ln(m^) 



^ ln(m^) - fff ln(rr&)S 



+ < 



942115 37 c 3403793 c 1 C i 643159 c _i_ 3758107 Jl , 10 ^2 i^/q^ , 7 
^18T - "g- ^3 - 19440 ^0 - 54 ^1 + ^2 + 116640 71 + "Tf ^ in Wj 

f§ ln(^) - ±f ln(^)S - ff (ln(t^)) 2 + 15 ln(w H ) \n(m 2 H ) + fj 
1 \n(m H )S + f (ln(m^ ^ 



I 61481 
995 



+ 1C3 

W H m H ) 



+ wh 



4522879 , 1329751 g , 120049 c i 39089 ^2 , 16039 w„„ \ 84613 i T1 / , -„ \o 
25920 + 38880 ^0 + ~W ^ 2 + ~2592" 71 + "96" [n { w H) ~ "H^Po 

™ ln«) - f ln«)S - 6 (ln(m| 



+ < 



+ 9 -§ (\n(w H )) 2 -§Hw H )\n(m 2 H 

239678009 , 15484933 g , 635063 
432000 129600 960 



74083 



35063 Q i 771181 Jl , 42518293 W„„ \ 9928 i„/„, \ c 

(ln(^)) 2 + 2f» ln(^) \n{m 2 H ) + ^ ln(m^) + 180 (ln(m^) 

16217189681 , 138228439 g , 1072891 c i 6308171 2 , 171948913 ] / \ 
2592000 48600 576 2 8640 ' 43200 1Ii l U7 ^/ 



604271 
648 



\n(w H )S + « (ln(^)) 2 + '-f Hwh) HO + ^ 



29 



D.2 Analytical results for the two-loop finite part of Z-boson 



A z - 



+ z H 



128 1 32 
2483 231 



_9_ 

32 



(\n(m 2 H 



m 



w 



519 
" 32 



ln(m^) + f| ln(m^ + f (HO) 



+ z 



H 



8311573 153 o 243 c 87651 q 47893 Jl , 459 /-(rt\ 399943 \„( „ \ 
~20736 8~ D3 T D ° 128~ D 2 ~~ "3456" 71 + ~16 ^ W 1728~ m \ Z H) 



+^ \n(z H )So + ^ (Hz H )) 2 + iffl ln(z H ) ln(mi) - = ln(m|) 



1595 



16 



+ 2 



86400 



10730119 3451 c 35739 c i 3167 _2 1173881 \„( „ \ , 517 l^/^ \o 
Q?T ^0 5n~ ^2 + TTisn 71 4320~~ ln { Z H) + ~ ^{Z H )O 

^ {Hz H )Y + ^ \n{z H )\n{m 2 H ) m , , , 



96 ' '" 80 2 1 1080 " 1320 
4 2 , 1743 , / \ L / 2\ 30383 



+ -2 



276774409 , 86473 o 3807 o _i_ 859 _2 168691 \ , 10323 \ o 



1296000 1 3200 



+ 



90121 
960 



(ln(z H )) 2 + ^ ln(z H ) ln(m H ) + ±|g3 + si 



+ z 



H 



13424129921 , 416057 c _ 452727 o , 1395311 2 _ 3016477 
J ii9n "2 ~r S n94n 10080 



14112000 



+ 



2400 



1120 



30240 



ln(z H ) 



7739 
40 



\n(z H )S + (ln(^)) 2 + Hz H ) \n(m H ) + 



6333137 
30240 



\n{m 2 H ) 



A? = z H 



2243 
' 288 



+ IlS + 2 -§S 2 + 2 ^n 2 + l \n(z H ) - \n(z H ) ln(m|) 



+§ ln(m|) - f ln(m^ 



+ Z H 

853 



4637683 , 589 c 12295 c _i_ aa a _i_ 157725 c i 99227 „2 , 11 2 l^/oN 
"^036^ + 12 ^3 - -34- ^0 + 00 ^1 + "64— ^2 + 71 + ~ 71 



-TC(3) 



864 



ln(^) - Ib(«h)5 - *f (ln(^)) 2 - *f In(z H ) ln(m|) 



203531 
864 



ln(m|) - ^ 1„«)S - ^ (ln(m^)) 



+ -2 



H 



4480681 _|_ 6445 a , 
43200 144 °° " t 



126999 
160 



5, 



1440 



7T^ + 



2160 



ln(^) - 2gZ \n(z H )S ( 



30 



144 



(Hz H ) f - 122 ln(z H ) \n(m 2 H ) + ^ \n(m 2 H ) - 99 ln(m 2 H )S - if (ln(m^ 



4 284392801 155881 c , 37971 c 661 _2 , 6123413 \ 3837 i„/, y \ c 

Z H 1296000 960~ ^° + ~W ^ 2 ~ ~18 71 + 43200 ln \ Z H> ~ ~ ^{Z H )^ 



11129 
96 



(ln(z H )) 2 - f ln(z H ) \n(m 2 H ) - |f ln(m^) - f (ln« 



^5 6259218463 1671301 g , 685701 c 5834099 Jl 270677 „ \ 25747 \n 
Z # 3969000 1800 D ° 224 ° 2 30240 77 6048 m l Z ^J 30 "H^PO 

-ffP (ln(^)) 2 + if ln(^) ln«) + « ln«) 



2ff 



3307 



, 7() £ So + 1 S 2 + 1 vr 2 + if ln(^) - g ln(^) ln(m|) - f ln(m^) 
+ 1 | ln(m|)5 - | (ln(m|)) S 



14576777 2923 q , 110455 c 1 7fi C 349773 g 111577 ^2 121 2 1t /q\ 
^0736 72" ^3 + 30 ~ I'D £1 T^T ^2 - ~J^T 71 ~ ~ 71 [n ( 6 ) 



+ C(3) + ff In(^) - g ln(* fl )So + I (M^)) 2 + W H*h) Hm 2 H ) 
+ iif|i ln(m|) - M l n «) 5o + |f (l n(m |) 



fZ 



452509 , 6145 



3200 



I b!45 c i 2913 c i 



2 913 
80 



2279 Jl 251521 



270 



7T 



1440 



ln(^) 



1943 
12 



ln(z H )S 



5623 
' 288 



(\n(z H )) 2 + f ln(z H ) ln(m^) - 2gi In(m^) + \n(m 2 H )S + f (ln(m^ 



106557211 3265133 
9600 



144000 



+ 



41949 
160 



S 2 + 



65291 Jl 16156451 



1440 



7T 



86400 



Mz H ) + ^i ln(^)S- 



+ ™ (ln(^)) 2 + ^ ln(^) Hm 2 H ) + ^ ln(m ^) + f ( ln(m |; 



160 

2 



109475704943 , 1 5528551 g , 1335741 o _i_ 9380087 2 , 34762379 
7200 1120 2 30240 ^ 151200 



25401600 



+ 



]n(z H ) 



188237 
120 



5529703 
30240 



ln(m 2 H ) 



+ 



77 g , 243 c _i_ 
144 ^0 + — ^2 + 



95 
108 



7T 2 + fi ln(^) - I \n(z H ) ln(m^) - if \n(m 2 H 



-ii Hm 2 H )S -\ (Hml))' 



31 



7 2 
Z H 



6479761 , 365 c 11525 c _i_ 457 c _i_ 102293 c _i_ 21419 Jl , 377 Ji n /q\ 

C(3) - 6 W H*h) + !f In(^)5 + I (Hz H )) 2 + f \n(z H ) ln(m^) 



7 3 



ln(mi) + ^ ln«)S + f (ln« 

315263 160771 c i 92279 c 3983 ^2 , 187433 „ \ 12461 i^/^ No 
13200" 2592" D ° + -430" D 2 - 4320 n + "Ueo - ln l^J 2l6~ m { z H)O 

f (ln(^)) 2 - f ln(^) ln(m^) + 2g° ln(rr&) - l -f ln(rr&)S - f (ln« 

275866331 3508117 q , 13353 c 1403 Jl , 1202389 i^/^ N 2932 \ 

1296000 10800" D + D 2 90~ ff + 10800 m \ Z H) 15- ^{ z H)O 

(ln(^)) 2 + ^ ln(^) ln(m|) + ^ ln(m|) + 9 (ln(m|)) 

178942516547 _ 871139 c i 8825401 c _ 6422629 2 _ 5656397 i / N 
63504000 324 °° ^ 3360 2 30240 50400 LLL \ Z H ) 

™ - ™ (ln(^)) 2 + W W**) ln«) + ™ ln«) 



7 4 



7 I 371 c I 1215 r. , 869 _2 , 133 1ti /„ \ 19 
9 + 864' b + ^2~^2 +8647r + ^{Z H ) - jg 



if ]n(z H ) H™ 2 h) + ^ H™>%) 



ln(m^)S - 



2 12699875 , 793 c i 1537817 c _ 365 o _ 22477 
H 20736 216 3 ' 2592 °° 9 1 128 

+ 3 -m C(3) + ff In(^) - H In(z fl )5 + | (ln(^)) 2 + W ln(** 
+ H ln«) - H ln«)S + ff (ln«)) ?1 

3 3647159 , 142763 o i 23303 c i 4529 2 _ 1921 
H 25920 7776 *~ 96 2592 " 144 



365 c _ 22477 c _ 1229275 
" ° l 128 2 31104 

r) ln(m^) 



3647159 , 142763 o , 23303 c _i_ 4529 J2 1921 1ti /„ \ , 743 i^/^ \o 
25920 + ^776~ ^° + "96" ^ 2 + 2592 71 ~ ~U4 ln K Z H) + 324 lr H^Po 

%g (\n(z H )) 2 - f \n(z H ) Hml) + *g ln(m^) + g ln(m|)5 - f (ln(m|)) 

696762899 1 36167453 q i 832139 c 1 38173 2 1862941 1 / \ , 56477 l / \o 
1296000 "t" 259200 °° 960 ° 2 8640 71 17280 ln \ z H ) "I" i 440 m l^^0 

+ ff (lnM) 2 + 92 lnNlnK) + ^f 

430433439517 , 343121543 o 4. 3744143 c 1 49949 2 _ 54474521 1 / \ 
127008000 ^ 194400 ^ 1120 2 *~ 672 151200 m v^-ff/ 



Z H 



2 -WMm 2 H ) + -, 



45 



(ln(m^)) 2 



z 5 



32 



2240561 Mz H )S + {Hzh )Y + f Hzh) H m 2 H ) + «i ln(m^) 



3240 "n^«7^u i 60 \m\~n/j i 24 "V'fl; *"V"HJ ' 1Q080 



fi+ 2 l^o + l^ + ifvr 2 + | ln(^) - | ln(^) lnK) 



+ ff ln(m^) - | Hm 2 H )S - § (ln(m 



+ ^ 



33773 1_ g 2421059 g _ 103 g , 219449 q , 1585397 2 , J_ 2 i /o\ 

384 108 3 19440 54 1 ^ 2880 2 233280 ~ l ~ 27 m v°/ 



-IC(3) 



72 

13429 
1440 



1440 



ln(zff) - if ln(^)S + (ln(^)) 2 + ^ ln(z H ) ln(m^ 



+ -2 



H 



+ 



32327203 , 90853 g , 202459 g 
259200 ~ t ~ 38880 ' 720 ■ 



17237 ^2 33671 



6480 



7T 



4320 



ln(^) - §g ln( Zff ) 5l 



-1 (ln(^)) 2 - f ln^)ln(m|) + iff ln^ 2 ,) - £ ln^)^ - 20 (ln(r< 



+ 4 



81023989 727891 



216000 43200 



S + 



164217 
160 



S 2 + 



10243 _2 78513 



1440 



7T — 



1600 



ln(^) - §§ ln(z H )S 



+ ff§ (ln(^)) 2 + f ln(^) ln^ 2 ,) + fp ln(m^) + f (lnKr 



+ -2 



23064502157 _ 107194699 
194400 



31752000 



S + 



19739897 
5040 



s 2 + 



361861 J2 , 4365419 



15120 



7T^ + 



151200 



In (z H ) 



« InW^o + H 1 (ln(^)) 2 + T In(^) ln(m*,) + ™ In(m^) 



References 

[1] D. Abbaneo et al. [ALEPH, DELPHI, L3 and OPAL Collaborations, LEP Electroweak 
Working Group, SLD Heavy Flavor and Electroweak Groups], |hep-ex/0 112021 . 

[2] M. Consoli, W. Hollik, F. Jegerlehner, in Z Physics at LEP1, 
eds. G. Altarelli et al., CERN 89-08 (1989) 1. 
For a more recent review see 

D. Bardin and G. Passarino, The Standard Model in the Making, 
Oxford, UK: Clarendon, 1999. 

[3] J. A. Aguilar-Saavedra et al. [ECFA/DESY LC Physics Working Group Collaboration], 
|hep-ph/01063T5] 



33 



[4] S. Eidelman and F. Jegerlehner, Z. Phys. C67 (1995) 585; 
M. Steinhauser, Phys.Lett. B429 (1998) 158. 

[5] J.J. van der Bij and M. Veltman, Nucl. Phys. B231 (1984) 205; 
J.J. van der Bij, Nucl. Phys. B248 (1984) 141; 
J.J. van der Bij and F. Hoogeven, Nucl. Phys. B283 (1987) 477; 

A. Djouadi, Nuovo Cim. A100 (1988) 357; 

M. Consoli, W. Hollik and F. Jegerlehner, Phys. Lett. B227 (1989) 167; 

B. A. Kniehl, Nucl. Phys. B347 (1990) 86; 

F. Halzen and B.A. Kniehl, Nucl. Phys. B353 (1991) 567; 
J. Fleischer et. al, Phys. Lett. B293 (1992) 437; 

K. G. Chetyrkin, A. Kwiatkowski and M. Steinhauser, Mod. Phys. Lett. A8 (1993) 
2785; 

G. Degrassi, Nucl. Phys. B407 (1993) 271; 

G. Buchalla and A. J. Buras, Nucl. Phys. B398 (1993) 285; 

R. Barbieri et al, Phys. Lett. B288 (1992) 95; B312 (1993) 511(E); 

Nucl. Phys. B409 (1993) 105; 

J. Fleischer, O.V. Tarasov and F. Jegerlehner, Phys. Lett. B319 (1993) 249; 
L. Avdeev et al., Phys.Lett. B336 (1994) 560; B349 (1995) 597(E); 
K. G. Chetyrkin, J.H. Kiihn and M. Steinhauser, Phys. Lett. B351 (1995) 331; 
K.G. Chetyrkin, J.H. Kiihn and A. Kwiatkowski, Phys. Lett. B282 (1992) 221; 
Phys. Rept. 277 (1996) 189; 

A. Czarnecki and J. H. Kiihn, Phys. Rev. Lett. 77 (1996) 395; 

R. Harlander, T. Seidensticker and M. Steinhauser, Phys. Lett. B426 (1998) 125; 

J. Fleischer et al., Phys. Lett. B459 (1999) 625; 

G. Degrassi, S. Fanchiotti and P. Gambino, Int. J. Mod. Phys.AlO (1995) 1377 

G. Degrassi et al., Phys. Lett. B350 (1995) 75; 

T. van Ritbergen and R.G. Stuart, Phys. Rev. Lett. 82 (1999) 488; 

D. Bardin et al., Comput. Phys. Commun. 133 (2001) 229. 

[6] G. Degrassi et al., Phys. Lett. B350 (1995) 75. 

[7] G. Degrassi, P. Gambino and A. Vicini, Phys. Lett. B383 (1996) 219; 
G. Degrassi, P. Gambino and A. Sirlin, Phys. Lett. B394 (1997) 188; 
G. Degrassi et al., Phys. Lett. B418 (1997) 209; 
S. Bauberger and G. Weiglein, Phys. Lett. B419 (1998) 333; 
A. Freitas et al., Phys. Lett. B495 (2000) 338. 

[8] G. Degrassi and P. Gambino, Nucl. Phys. B567 (2000) 3. 

[9] T. van Ritbergen and R.G. Stuart, Phys. Rev. Lett. 82 (1999) 488; 
Nucl. Phys. B564 (2000) 343; 

M. Steinhauser and T. Seidensticker, Phys. Lett. B467 (1999) 271. 
[10] F. Jegerlehner, Nucl. Phys. Proc. Suppl. 37B (1994) 129. 



34 



[11] A. Sirlin, in Proc. of the 19th Intl. Symp. on Photon and Lepton Interactions at High 
Energy LP99 ed. J.A. Jaros and M.E. Peskin, Int. J. Mod. Phys. A 15S1 (2000) 398. 

[12] V.A. Smirnov, Commun. Math. Phys. 134 (1990), 109; Mod. Phys. Lett. A10 (1995) 
1485; Renormalization and asymptotic expansions (Bikrhauser, Basel, 1991). 

[13] P. Gambino, A. Sirlin and G. Weiglein, J.High Energy Phys. 04 (1999) 025. 

[14] G. 't Hooft and M. Veltman, Nucl. Phys. B50 (1972) 318; 
D.A. Ross and J.C. Taylor, Nucl. Phys. B51 (1973) 125; 

D.Yu. Bardin, P. Khristova and O.M. Fedorenko, Nucl. Phys. B175 (1980) 435; 

Nucl. Phys. B197 (1982) 1; 

S. Sakakibara, Phys. Rev. D24 (1981) 1149; 

K. Aoki et al., Prog. Theor. Phys. Suppl. 73 (1982) 1; 

M. Bohm, H. Spiesberger and W. Hollik, Fortsch. Phys. 34 (1986) 687; 

A. Denner, Fortsch. Phys. 41 (1993) 307. 

[15] A. Sirlin, Phys. Rev. D22 (1980) 971. 

[16] J. Fleischer and F. Jegerlehner, Phys. Rev. D23 (1981) 2001; 
F. Jegerlehner, in "Radiative Corrections in SU(2)l x £7(1)", 
eds. B.W. Lynn, J.F. Wheater, World Scientific Publ., Singapore, 1984. 

[17] J. Fleischer, O. V. Tarasov and M. Tentyukov, Nucl. Phys. Proc. Suppl. 89 (2000) 112. 

[18] C. Ford, I. Jack and D.R.T. Jones, Nucl. Phys. B387 (1992) 373; B504 (1997) 551(E); 
A.I. Davydychev and J.B. Tausk, Nucl. Phys. B397 (1993) 123. 

[19] G. Weiglein, R. Scharf and M. B6hm, Nucl. Phys. B416 (1994) 606. 

[20] D.J. Broadhurst, J. Fleischer and O.V. Tarasov, Z. Phys. C60 (1993) 287; 
F.A. Berends et. al., Z. Phys. C63 (1994) 227; 

S. Bauberger, F.A. Berends, M. Bohm and M. Buza, Nucl. Phys. B434 (1995) 383. 

[21] R. Scharf and J.B. Tausk, Nucl. Phys. B412 (1994) 523. 

[22] F.A. Berends and J.B. Tausk, Nucl. Phys. B421 (1994) 456; 
A. Ghinculov and J.J. van der Bij, Nucl. Phys. B436 (1995) 30; 
S. Bauberger and M. Bohm, Nucl. Phys. B445 (1995) 25. 

[23] L. Baulieu and R. Coquereaux, Ann. Phys. 140 (1982) 163. 

[24] P. Nogueira, J. Comput. Phys. 105 (1993) 279. 

[25] J. Fleischer and M.N. Tentyukov, Comp. Phys. Commun. 132 (2000) 124. 

[26] J. A. M. Vermaseren, Symbolic manipulation with FORM, 
Amsterdam, Computer Algebra Nederland, 1991. 



35 



[27] O.V. Tarasov, Nucl. Phys. B502 (1997) 455; Phys. Rev. D54 (1996) 6479. 

[28] J. Fleischer, M.Yu. Kalmykov and O.L. Veretin, Phys. Lett. B427 (1998) 141. 

[29] L.V. Avdeev et al, Nucl. Ins. Meth. A389 (1997) 343; 
Comp. Phys. Commun. 107 (1997) 155. 

[30] J. Fleischer et al., Phys. Lett. B459 (1999) 625. 

[31] M. Veltman, Diagramatica, Cambridge University Press, 1994. 

[32] J. Fleischer and O.V. Tarasov, Comp. Phys. Commun. 71 (1992) 193; 

J. Fleischer and M. Yu. Kalmykov, Comp. Phys. Commun. 128 (2000) 531. 

[33] J. Fleischer et al, Nucl. Phys. B539 (1999) 671; B571 (2000) 511(E). 

[34] J. Fleischer, A.V. Kotikov and O.L. Veretin, Nucl. Phys. B547 (1999) 343; 
J. Fleischer, M.Yu. Kalmykov and A.V. Kotikov, Phys. Lett. B462 (1999) 169; 
B467 (1999) 310(E); 



A.I. Davydychev and M.Yu. Kalmykov, Nucl. Phys. B605 (2001) 266 (|hep-th/0012189D 



[35] R. Tarrach, Nucl. Phys. B183 (1981) 384; 

A.S. Kronfeld, Phys. Rev. D58 (1998) 051501. 

[36] P. Gambino and PA. Grassi, Phys. Rev. D62 (2000) 076002. 

[37] G. 't Hooft and M. Veltman, Nucl. Phys. B44 (1972) 189; 
C.G. Bollini and J.J. Giambiagi, Nuovo Cim. B12 (1972) 20; 
J.F. Ashmore, Lett. Nuovo Cim. 4 (1972) 289; 
G.M. Cicuta and E. Montaldi, Lett. Nuovo Cim. 4 (1972) 329. 

[38] D.J. Gross and F. Wilczek, Phys. Rev. D8 (1973) 3633. 

[39] G. 't Hooft, Nucl. Phys. B61 (1973) 455. 

[40] D.R.T. Jones, Nucl. Phys. B87 (1975) 127; Phys. Rev. D25 (1982) 581; 
M.E. Machacek and M.T. Vaughn, Nucl.Phys. B222 (1983) 83. 

[41] D.J. Broadhurst, Z. Phys. C54 (1992) 599. 

[42] N. Gray et al., Z. Phys. C48 (1990) 673; 

L.V. Avdeev and M.Yu. Kalmykov, Nucl. Phys. B502 (1997) 419. 

[43] M. Veltman, Acta Phys. Polon. B8 (1977) 475. 

[44] G. 't Hooft and M. Veltman, Nucl. Phys. B153 (1979) 365. 

[45] U. Nierste, D. Muller and M. B6hm, Z. Phys. C57 (1993) 605. 



36 



[46] A.I. Davydychev, Proc. Workshop "AIHENP-99", Heraklion, Greece, April 1999 
(Parisianou S.A., Athens, 2000), p. 219 ( |hep-th/9908032|) ; Phys. Rev. D61 (2000) 
087701. 

[47] A.I. Davydychev and M.Yu. Kalmykov, Nucl. Phys. B (Proc. Suppl.) 89 (2000) 283 
(|hep-th/0005287| ); Nucl. Phys. B605 (2001) 266 flhep-th/0012189|) . 

[48] L. Lewin, Polylogarithms and associated functions (North-Holland, Amsterdam, 1981). 

[49] F. A. Berends, A. I. Davydychev and V. A. Smirnov, Nucl. Phys. B478 (1996) 59. 

[50] E.E. Boos and A.I. Davydychev, Theor. Math. Phys. 89 (1991) 1052. 

[51] D.J. Broadhurst, |hep-th/9604128| ; Eur. Phys. J. C8 (1999) 311; 
J. Fleischer and M. Yu. Kalmykov, Phys. Lett. B470 (1999) 168; 
M.Yu. Kalmykov and O. Veretin, Phys. Lett. B483 (2000) 315. 

[52] A. Freitas, W. Hollik, W. Walter, G. Weiglein, "Electroweak two-loop corrections to 
the Mw — Mz mass correlation in the Standard Model" , [hep-ph/0202131| . 



37 



-- 1 - 

.... 2 - 


loop 
loop 


/ 


sum 




/ 

/ 

/ / 
/ / 
/ / 

/ / 
/ / 
/ / 
/ / 

/ jS 






/ - - - . 
/ " - - . . 






/ 

s ■ — 






^ - 

^^^^ 


100 


150 


200 250 300 



M H (GeV) 

Figure 3: One- and two-loop corrections to the relation A z = m|(M z )/M| — 1 as a function 
of the Higgs mass M H for intermediate Higgs masses. 
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Figure 4: One- and two-loop corrections to the relation A z = m|(M z )/Mf — 1 as a function 
of the Higgs mass M H for heavy Higgs masses. 
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Figure 5: One- and two-loop corrections to the relation Ay/ = vr^iM^jM^ — 1 as a 
function of the Higgs mass m H for intermediate Higgs masses. 
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Figure 6: One- and two-loop corrections to the relation A w = rri^ v (M w ) / — 1 as a 
function of the Higgs mass m H for heavy Higgs masses. 
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Figure 7: The dependence on the number of coefficients of the expansion ( |5.4|) used for the 
evaluation of the two- loop corrections. We show 5 = — (p-z + fig fig') /M§ (see |5.3|) as 
a function of the Higgs mass. The dotted, dashed, dot-dashed and full lines show results 
obtained with the first one, two, three and all calculated (six) coefficients, respectively. 



0.002- 

0.001- 

0.000- 
^-0.001- 
-0.002- 
-0.003- 
-0.004- 

100 150 200 250 300 

M H (GeV) 

Figure 8: The dependence on the number of coefficients of the expansion with respect to 
sin 2 8w of the two-loop corrections. We show 5 = — (fl^ + ) l^z ( see j[lD as 

a function of the Higgs mass. The dotted, dashed, dot-dashed and full lines show results 
obtained with the first one, two, three and all calculated (six) coefficients, respectively. 
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Figure 9: One- and two-loop corrections to A s = sin 2 9^p/ sin 2 9yy — 1 (see ( |5.8|) ) as a 
function of the Higgs mass m# for intermediate Higgs masses (/i = Mz)- 
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Figure 10: One- and two-loop corrections to A s = sin 2 6yp/ sin 2 9^ — 1 (see (|5.8|)) as a 
function of the Higgs mass mn for heavy Higgs masses (/x = Mg). 
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